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ABSTRACT 

We report the first attempt to observe the secondary eclipse of a transiting extra- 
solar planet at radio wavelengths. We observed HD 189733 b with the Robert C. 
Byrd Green Bank Telescope of the NRAO over about 5.5 hours before, during and 
after secondary eclipse, at frequencies of 307 - 347 MHz. In this frequency range, we 
determine the 3-f7 upper limit to the flux density to be 81 mJy. The data are consistent 
with no eclipse or a marginal reduction in flux at the time of secondary eclipse in all 
subsets of our bandwidth; the strongest signal is an apparent eclipse at the 2-a level 
in the 335.2 - 339.3 MHz region. Our observed upper limit is close to theoretical 
predictions of the flux density of cyclotron-maser radiation from the planet. 

Key words: planetary systems - stars: HD 189733 - radio continuum: stars - stars: 
magnetic fields - masers 
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1 INTRODUCTION 

In the solar system, all the magnetised planets are known to 
emit at radio wavelengths. One of the principal mechanisms 
producing this is electron-cyclotron maser radiation, which 
produces emission at the local gyrofrequency, /g, given by 



/g(MHz) = 2.8B(G), 



(1) 



where B is the planetary magnetic field strength. 

Of the solar system planets, Jupiter has the strongest 
radio emission, outshining the (quiescent) Sun by several or- 
ders of magnitude at frequencies of about 30 MHz. Jupiter 
was also the first planet discover ed to emit at radi o wave - 
lengths, detected at 22 MHz bv iBurke fc FrankUnI (|l95d ). 
since the other solar system planets (including Earth) emit 
at frequencie s below the Earth's io nospheric cut-off at ~ 
5-10 MHz l|Griefimeier et al.ll2005h . 

It has long been suggested that analogous radio emis- 
sion might arise in an extra-solar planetary system, and that 
it may be possible to detect this radiation from Earth. In- 
deed, several attempts have been made to detect radio emis- 
sion from an extra-solar planetary system, beginning with 



* E-mail: amssOst-and. ac.uk 

t Scottish Universities Physics Alliance 



observations in the 197 0s and 80s before any extra-solar 
planets were known (see iGriefimeier et al. I I2OO6I for a list of 
observational campaigns) . More recently, with the discovery 
of hundreds of extra-s olar planets starting with 51 Peg b 
(|Mavor &: Quelo3ll995l ). several further attempts have been 
made to detect radio emission from these planets, none suc- 
cessfully to date. 

Observations of the radio emission from Jupiter have 
allowed the planet's magnetic field strength to be inferred 
from the high-frequency cut-off (around 40 MHz). The emis- 
sion with this highest frequency is generated in the region 
with the highest magnetic field strength (equation [ij , close 
to the planet. Observations of the radio spectrum of hot 
Jupiters could provide similar information about the plane- 
tary magnetic field, composition and rotation. 

Although Jupiter's emission is several orders of mag- 
nitude stronger than the (quiescent) Sun's at certain fre- 
quencies, the analogous emission from extra-solar systems is 
expected to be very faint because of the distance involved. 
Consequently, the planetary systems most favoured by ob- 
servers tend to be those within a few pc of Earth. 

A new theoretical model (jjardine fc Collier CameronI 

l2008h predicts the radio flux from a range of exoplanets, both 
those in large orbits that are embedded in the stellar wind, 
and those in close orbits that lie within the stellar magne- 
tosphere. Reconnection between the stellar and planetary 
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magnetic fields provides the pool of accelerated electrons 
necessary for the electron-cyclotron maser instability. This 
model reproduces the observed emission from solar-system 
planets and predicts that the radio flux from exoplanets in 
close orbits should scale as 



Si, oc 



N 



B. 



1/3 



(2) 



where A'^ is the coronal density and is the surface field 
strength of the parent star. 

In particular, this model predicts the radio flux for 
HP 189733 b, a h ot- Jupiter planet orbiting a K2 dwarf 
Bouchv et"ai]2005lV For the observed held strength of 40 G 
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200?! ') and a solar coronal density, this gives 



15mJy, at a frequency determined by the planetary mag- 
netic held strength. It is likely, however, that the coronal 
density will be greater than the solar value, since HD189733 
rotates more rapidly than the Sun . If we assume scalings 
in the range N oc n°-^ to N cx ft (|Unruh fc Jardind [l997l . 
Ilvanova fc TaarnI |20 03l) then for the rotation rate of 11.73 
days appropriate for HD189733, the flux rises to 39 - 72 
mjy respectively. 

In this work, we make use of the fact that HD 189733 b 
transits its host star to observe the system during secondary 
eclipse, when the planet passes behind the star. Detecting 
the secondary eclipse at radio wavelengths would provide un- 
ambiguous evidence that the emission is of planetary, rather 
than stellar, origin. These observations (described in Section 
[2|, conducted on a bandwidth centred at 327 MHz, are used 
to determine new upper limits to the radio flux from HD 
189733 b (Section This relatively low frequency was cho- 
sen because the planet could emit at this frequency given 
a plausibly strong magnetic held. The signiflcance of these 
upper limits and the prospects of future work are discussed 
in Section |4l 



2 OBSERVATIONS 

We observed the HD 189733 system fo r ~5.7 hours, encom- 
passing a 1.827 hour (|Winn et al.ll2007l ) secondary eclipse on 
2007 April 21, when the star was at altitudes between 45° 
and 74°. The observations were conducted using the Spec- 
tral Processor at PFl of the Robert C. Byrd Green Bank 
Telescope (GBT) of the National Radio Astronomy Obser- 
vatory (NRAO) 

The dual polarisation observations span the frequency 
range 307 MHz to 347 MHz and consist of a series of spectra 
comprising 1024 frequency channels, each with an integra- 
tion time of one second. Spectra were grouped into scans, of 
duration 120 s, and scans of the target (ON) were alternated 
with off-target scans (OFF) of a patch of sky whose centre 
is 2° South of the target in declination. 

The reason for choosing to make spectral-line observa- 
tions, despite the broad-band nature of the expected signal, 
is to mitigate against radio frequency interference (RFI). 
Much RFI is narrow-band in nature, so observing with a 
large number of frequency channels allows channels affected 
by RFI to be excluded (see Section [5T|> . 

Because of the large beam size at this frequency, we 
expect there to be other radio sources within the beam. The 



WENSS 300 MHz survey has a sensitivity of 18 mJy, but 
does not extend below +30° in declination and so does not 
cover our target star; it can, however, give us an idea of how 
many objects we can expect in a typical GBT beam (FWHM 
= 36'). We looked in the WENSS catalogue at ten beam 
pointings spaced between 20 and 21 hours in RA, at +31° 
in declination. We flnd an average of 4.8 objects within 20' 
of the beam centre; these objects have a median integrated 
flux of 121 m Jy. The Texas Survey of discrete radio sources 
at 365 MHz (|Douglas et al.|[l996l ) is less sensitive (90 per 
cent complete to 400 mJy and 80 per cent to 250 mJy), 
but it does cover HD 189733. The Texas survey includes an 
object with a flux density of 550 mJy 18.8' from the centre 
of our ON beam, and two objects of flux density 581 mJy 
and 241 mJy, respectively 3.6' and 16.1' from the centre of 
our OFF beam. 

Given that we are searching for a source exhibiting a 
speciflc, known variabilitjQ, the other radio sources in our 
beam are unlikely to be a problem. This is because in order 
to mimic the eclipse we are looking for, any variable object 
would not only have to vary, but would have to produce 
a large (~ 10 per cent or greater) eclipse-like signal on a 
speciflc, relatively short, timescale. 

Each observation was calibrated using a noise tube cal- 
ibrator, but we do not correct for atmospheric opacity, since 
such effects are negligible in the 1 m wavelength regime. 
Given the beam size and pointing stability of the GBT, we 
do not need to correct for pointing errors. An increase in raw 
flux was observed in both our target and off-target beams at 
around 2454214.0 JD which corresponds to 8 am local time 
and is probably caused by an increase in radio-frequency 
interference. 



3 ANALYSIS 

3.1 Production of lightcurves 

Initially, the data were visually inspected by stacking a num- 
ber of spectra into a time series. A typical set of stacked 
spectra are shown in Figure [TJ using colour to indicate flux. 
The most noticeable features in these stacked spectra are 
(i) the sharp 'spikes' spanning just one or two frequency 
channels and (ii) an oscillation with a period of around 20 s 
across all frequencies. Since both these effects are observed 
in both the ON and OFF scans, we conclude that they are 
due to radio-frequency interference. 

In order to determine which frequency channels suffer 
from intermittent 'spiky' RFI, the variance of each of the 
1024 channels over the whole dataset was calculated (Fig- 
ure 



2(a) I. After noticing a small number of pairs of scans 



suffer from a high level of RFI over a broad range of fre- 
quencies, we identifled and excluded a total of seven pairs of 
scans from further analysis. It is believed that the RFI in at 
least some of these scans was caused by electrical interfer- 
ence in the control room of the GBT. The variance in each 
frequency channel after the exclusion of these scans (Figure 



^ Using the period and transit duration of IWinn et al.1 ll2007t) . 
and the time of mid-eclipse from lDeming et al. I boost) , and asso- 
ciated uncertainites, we conclude that the eclipse time is known 
to within 200 s. 
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(a) 



Figure 1. Series of time-stacked spectra for a typical pair of scans 
(ON, lower half of figure; OFF, upper half), each of duration 
120 s. Time increases from bottom to top; frequency from left 
(307 MHz) to right (347 MHz). Colours are indicative of raw flux 
counts, red (light) represents high flux and blue (dark) low flux. 



2(b) I shows only narrow band spikes of RFI, which is easily 
identified and removed. Note that the same analysis of the 
raw OFF scans reveals spikes of high flux density in identi- 
cal frequency channels, confirming that RFI, and not fiaring 
of the target star, is responsible. A total of 177 of the 1024 
channels were in this way identified as containing RFI and 
excluded, resulting in a variance plot without spikes (Figure 

Ml- 

In order to determine which parts of the remaining 
bandpass to use to construct lightcurves and search for an 
eclipse signal, we calculated the standard deviation of each 
frequency channel of the raw ON lightcurve, which is nor- 
malised by dividing by the mean raw flux (flgure [5]) . With 
the aid of this figure, four regions of the bandpass that show 
less variation than surrounding regions were identified. We 
number these four subsets of the bandpass (i) - (iv) and 
construct lightcurves from the frequency channels in each of 
these bandpasses, as well as one from all of the remaining 
847 frequency channels. 

In each case, single, uncalibrated ON and OFF 
lightcurves were created by taking the mean of the rele- 
vant channels at 1 s intervals. The periodic RFI identified 
in Figure [T] earlier is clearly visible in plots of these ON and 
OFF lightcurves spanning only a few minutes (Figure |3]). A 
rectified sine wave was used to fit this periodic modulation, 
which we believe to be caused by a distant radar system. 
The ON and OFF lightcurves for the whole bandpass, after 
the rectified sine wave has been subtracted, are shown in 
figure |4l 

After subtracting the rectified sine wave from each scan 
(ON and OFF), the lightcurves were median-binned on the 
duration of a scan, and the ON scan was converted to a flux 
density, S^, in Jy, according to 



1 

r 



{ON - OFF) 
OFF 



(3) 



where ON and OFF are the raw fluxes of the ON and OFF 
scans respectively, Tsys is the system temperature, and F — 
2.0 KJy~^ is the telescope sensitivity. The resulting binned, 
calibrated lightcurves are shown in Figure [6] Although the 
residuals in our fit to the periodic RFI result in errors on 
individual points of magnitude ~ 100 mjy, this should not 
affect our ability to detect a signal, given our number of data 
points. 
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Figure 2. Variance in the raw flux for the ON scans taken on 2007 
April 23. The top panel includes all 66 ON scans and all frequency 
channels. The middle panel is the same as the top panel, but 7 
noisy scans have been completely excluded. The bottom panel 
is the same as the middle panel, but with 177 noisy frequency 
channels removed. The four regions indicated in each panel are 
those parts of the bandwidth for which we generate lightcurves 
(see text). 



3.2 Measuring the eclipse depth 



We use the period (P — 2.2185733 days), epoch of transit 
{to = 2453988.80336 (HJD)), and transit duration {tiv - 
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Figure 3. A typical ON scan of duration 120 s sliowing tlie raw 
counts (averaged over all frequency channels) as a function of 
time (+). Also plotted are the rectified sine wave fitted to this 
scan (solid line) and the residuals of this fit (x). 
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Figure 5. Standard deviation of each channel of raw ON data, 
normalised by the mean flux of each channel as a function of fre- 
quency. The four regions indicated are those parts of the band- 
width for which we generate lightcurves (see text) . 
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Figure 4. Uncalibrated lightcurves ON (upper curve; denoted 
by +) and OFF (lower; x) target, after subtraction of periodic 
(rectifled sine wave) RFI. 



upper limit to the flux density for each of the five lightcurves 
of figure [6] are displayed in table [1] 

Over the whole bandwidth, we do not detect a signifi- 
cant eclipse, but we are able to place a 3-a upper limit of 81 
mjy to the flux density from the planet. Of the four selected 
subsets of the bandpass, we find insignificant brightening 
(consistent with a constant flux model) at the predicted time 
of secondary eclipse in regions (i) and (ii) . Eclipses with sig- 
nificances of 1.5-a and 2.0-a are found in regions (iii) and 
(iv), respectively. 

It is possible that these 'detections' of the radio eclipse 
are caused by scatter in the lightcurve, and are not real; 
we do not understand why it should be any easier to de- 
tect the eclipse signal in these particular narrow frequency 
ranges than in any other. We are unable to conclude with 
any confidence that we have detected the eclipse, although 
these offer the possibility that we are tantalisingly close to 
detecting the first radio emission from an extra-solar planet. 



ti = 1.827 hours) determined by I Winn et all l|2007l ) and the 
fact that the orbit of HD 189733 b is circular to identify the 
'low' points, £, observ ed during secondary ec lipse. We then 
follow the method of ICollier Cameron et ali (|2[)065 to cal- 
culate the depth of transit and the associated uncertainty. 
The weighted mean of the lightcurve is subtracted from all 
lightcurve points, Xi, and the inverse- variance weights, Wi, 
are scaled to give = 1 for ^ constant-flux model fitted to 
the out-of-eclipse data. The transit depth, 5, and its vari- 
ance, Var(5) are then calculated according to. 



and 



S = 



Var(5) 



(4) 



(5) 



respectively. With the null hypothesis that there is no 
eclipse, the n-cr upper limit to the flux density is given by 
nVVar(5). 

The depth of transit, associated uncertainty and the 3-cr 



4 DISCUSSION 

4.1 Interpretation of our results 

There are three aspects of planetary cyclotron maser emis- 
sion that might lead to a non-detection: (i) the planet does 
not radiate at the observed frequencies, (ii) any emission is 
insufficiently strong for us to detect and (iii) the emission is 
beamed out of the line of sight 

We consider first the case where the magnetic field 
strength of the planet does not allow radiation to be gen- 
erated at the observed frequency range. It seems unlikely 
that the magnetic field is too strong, since it is expected 
that the low-frequency tail of the spectrum will be fairly 
flat (|Griefimeier et all 120051 ). The fact that our higher fre- 
quency observations produce more encouraging results, how- 
ever, hints at the possibility that the magnetic field is too 
strong to produce emission at the low-frequency end of the 
observed bandwidth. 

Jupiter's cyclotron maser emission is observed from 
about 0.3 to 40 MHz, corresponding to a highest local value 
of the magnetic field strength of about 14 G. In order for us 



Radio observations of HD 189733 h 

Table 1. Measured eclipse depths and upper limits to the radio flux density from HD 189733 b in several parts of the bandwidth. 
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Bandpass 


Frequency 


Eclipse 


Standard deviation 


Significance 


3-(T upper limit to 




range (MHz) 


depth (mjy) 


of depth (mJy) 




flux density (mJy) 


all 


307.0 - 347.0 


-21.6 


27.1 


-0.8 cr 


81 


(i) 


308.6 - 312.7 


+15.6 


31.4 


+0.5 a 


94 


(ii) 


319.3 - 325.6 


+6.8 


26.5 


+0.3 a 


80 


(iii) 


329.8 - 332.7 


-50.3 


33.9 


-1.5 a 


102 


(iv) 


335.2 - 339.3 


-56.5 


28.0 


-2.0 a 


84 



(b) 



(d) 




Figure 6. Calibrated lightcurves of HD 189733 in different frequency ranges, (a) the whole bandpass, and (b), (c), (d) and (e), the four 
bandpasses (i), (ii), (iii) and (iv) respectively (see figure [Sj- Also plotted in each panel is the best fit to the data (sec text for details). 



6 A. M. S. Smith et al. 



to detect emission at 307 - 347 MHz, the emission would have 
to be generated in a region where the local field strength is 
~110 G. Since it is probable that Jupiter's emission orig- 
inates well above the cloud deck, it is possible that emis- 
sion from an extra-solar planet could be produced closer to 
the planetary surface, where the local field strength would 
be greater. However, if no emission is generated in regions 
where the magnetic field strength is about 110 G or stronger, 
all of the emission would occur at frequencies below our ob- 
serving bandwidth and no signal would be detected. 

The second case is that in which the planet may produce 
emission in the observed frequency range, but at an insuffi- 
cient intensity to be detected b y our observations. Accord- 
ing to the theoretical model of Ijardine fc Collier CameronI 
l|2008l ). a stellar coronal density between 2.0 and 2.5 times 
greater than the solar value implies the existence of a popu- 
lation of electrons sufficient to generate radiation that would 
be detected at our 3-cr confidence level. Therefore, the stel- 
lar coronal density must be at least twice solar in order for 
us to have been able to detect the emission, even if our ob- 
servations were conducted at the optimum frequency. 

Third, the emission may be beamed in a cone, in a man- 
ner similar to Jupiter's lo-controUed radiation. If this is the 
case it may be that the beam does not cross our line of 
sight and we would not observe any radio emission. Predict- 
ing this possible beaming effect is beyond the scope of our 
current model, however. 

Additionally, there exist potential observational reasons 
for the lack of detection. One possibility is that the integra- 
tion time of individual spectra was too long, at 1 s, given 
that Jupiter's emission displays variability on very short 
timescales. The emission from Jupiter that varies on the 
shortest (millisecond) timescales, however, is that controlled 
by the moon lo, whereas the emission not associated with l o 
varies on timescales of order Is or greater (jCarr et al.| | l983j ). 
Furthermore, even if there is short period variability of the 
emission from HD 189733b, we should still be able to ob- 
serve the mean level of this flux drop as the planet goes into 
eclipse. 

Another possibility is that the eclipse was masked by 
the variability of another radio source in either the ON or 
the OFF beam. This is discussed in Section^ where we sug- 
gest that it is unlikely that the eclipse signal is mimicked or 
masked by the variability of other objects, given the nature 
and timing of the expected signal. 



4.2 Comparison with previous observations 

The upper limits calculated in this work compare favourably 
with those de r ived in previous work on different systems. 
iBastian et all (j2000) ) observed seven extra-solar planetary 
systems with the Very Large Array (VLA), with typical 1 
(J sensitivities of around 50 mjy at 74 MHz an d 1-10 mjy 
at 333 MHz. More recently, iLazio et aP ()2004h observed 5 
systems with the VLA and derived 2.5 a upper limits of 218 
to 325 mJy. A further limit of 150 - 300 mJy using the VLA 
was placed on the radio fl ux from the extra-solar planet r 
Boo iLazio fc Farrellll2007^ . Tighter (2.5-a) upper limits of 
7.9 and 15.5 mJy were placed o n the emission from e Er i b 
and HD 128311b respectively bv lGeorge fc Stevens! (120071 ) at 
150 MHz using the Giant Metrewave Radio Telescope. Ours 



are the flrst eclipse observations, and they produce upper 
limits consistent with previous observations. 



5 CONCLUSIONS AND FUTURE PROSPECTS 

We have determined new upper limits to the radio flux den- 
sity from the transiting planet HD 189733 b. Our 3-a limit 
to the flux density over 307 - 347 MHz is 81 mJy. Of the four 
subsets of this bandpass selected for their noise properties, 
we detect no eclipse in two of them, but we detect marginal 
eclipses in the other two, the strongest of which is a 2-a de- 
tection. Further observations at the time of secondary eclipse 
in the future may enable us to improve the signal-to-noise 
further in order to determine whether the signal apparently 
detected is real or not. 

Ultimately, it may be that a more sensitive, low fre- 
quency instrument such as the forthcoming Low Frequency 
Array (LOFAR) is required to make the first detection of 
radio emission from the magnetosphere of an extra-solar 
planet. Not only will LOFAR be more sensitive than exist- 
ing instruments, but it will be able t o observe at frequencies 
of 10 - 240 MHz (|Farrell et al.ll2004l '). enabling the detection 
of emission from planets with weaker magnetic fields than 
at present. 



6 ACKNOWLEDGEMENTS 

AMSS acknowledges the financial support of a UK PPARC / 
STFC studentship. This research has made use of the VizieR 
catalogue access tool, CDS, Strasbourg, France. The authors 
wish to thank the anonymous referee, whose helpful com- 
ments led to improvements in the manuscript; and Keith 
Horne, for useful discussions. 



REFERENCES 

Bastian T. S., Dulk G. A., Leblanc Y., 2000, ApJ, 545, 
1058 

Bouchy F., Udry S., Mayor M., Moutou C, Pont F., Irib- 
arne N., da Silva R., Ilovaisky S., Queloz D., Santos N. C, 
Segransan D., Zucker S., 2005, AfcA, 444, L15 

Burke B. F., Franklin K. L., 1955, J. Geophys. Res., 60, 
213 

Carr T. D., Desch M. D., Alexander J. K., 1983, Phe- 
nomenology of magnetospheric radio emissions. Physics 
of the Jovian Magnetosphere, pp 226-284 

Collier Cameron A., PoUacco D., Street R. A., Lister T. A., 
West R. G., Wilson D. M., Pont F., Christian D., Clarkson 
W. I., Evans A., HasweU C. A., Helher C, Hodgkin S. T., 
Horne K., Irwin J., Kane S., Norton A., Skillen I., Wheat- 
ley P., 2006, MNRAS, submitted, astro-ph/0609418 

Deming D., Harrington J., Seager S., Richardson L. J., 
2006, ApJ, 644, 560 

Douglas J. N., Bash F. N., Bozyan F. A., Torrence G. W., 
Wolfe C, 1996, AJ, 111, 1945 

FarreU W. M., Lazio T. J. W., Zarka P., Bastian T. J., 
Desch M. D., Ryabov B. P., 2004, Planetary and Space 
Science, 52, 1469 

George S. J., Stevens I. R., 2007, MNRAS, 382, 455 



Radio observations of HD 189733 b 7 

Griefimeier J. M., Motschmann U., Glassmcicr K. H., Mann 
C, Rucker H. O., 2006, in Arnold L., Bouchy F., Moutou 
C, cds, Tenth Anniversary of 51 Peg-b: Status of and 
prospects for hot Jupiter studies The potential of exo- 
planetary radio emissions as an observation method, pp 
259-266 

Griefimeier J.-M., Motschmann U., Mann G., Rucker H. O., 

2005, A&A, 437, 717 
Ivanova N., Taam R. E., 2003, ApJ, 599, 516 
Jardinc M., Collier Cameron A., 2008, A&A, in press, 

arXiv:0808.2889 
Lazio T. J. W., Farrell W. M., 2007, ApJ, 668, 1182 
Lazio T. J. W., FarrcU W. M., Dictrick J., Grccnlccs E., 

Hogan E., Jones C, Hcnnig L. A., 2004, ApJ, 612, 511 
Mayor M., Qucloz D., 1995, Nat, 378, 355 
Moutou C, Donati J.-F., Savalle R., Hussain G., Alecian 

E., Bouchy F., Catala C, Collier Cameron A., Udry S., 

Vidal-Madjar A., 2007, A&A, 473, 651 
Unruh Y. C, Jardine M., 1997, A&A, 321, 177 
Winn J. N., Holman M. J., Henry G. W., Roussanova A., 

Enya K., Yoshii Y., Shporer A., Mazeh T., Johnson J. A., 

Narita N., Suto Y., 2007, AJ, 133, 1828 

This paper has been typeset from a file prepared 

by the author. 



